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Ab initio calculations using the local spin density approximation and also including the Hubbard 
U have been performed for three low energy configurations of the interface between LaAlOa and 
Ti02-anatase. Two types of interfaces have been considered: LaO/Ti02 and A102/TiO, the latter 
with Ti-termination and therefore a missing oxygen. A slab-geometry calculation was carried out 
and all the atoms were allowed to relax in the direction normal to the interface. 

In all the cases considered, the interfacial Ti atom acquires a local magnetic moment and its 
formal valence is less than +4. When there are oxygen vacancies, this valence decreases abruptly 
inside the anatase slab while in the LaO/Ti02 interface the changes are more gradual. 
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I. INTRODUCTION 

Complex oxide heterostructures have been the subject 
of many recent papers, both experimental and theoret- 
ical, as their interesting interface properties promise to 
pave the way towards novel electronic devices. With the 
current available experimental techniques oxide thin films 
can be produced with a high degree of cristallinity and 
the electronic structure of their surfaces and interfaces 
can be precisely determined. The most studied system 
of this type has been LaAlOa/SrTiOa (LAO/STO), that 
consists of two perovskite structures stacked along the 
(001) direction^, in which they present alternate layers of 
LaO, A102, Ti02 and SrO with a very small lattice mis- 
match between them. A similar but less studied interface 
is that of LAO with Ti02-anatase, that presents an even 
smaller lattice mismatch, namely less than 0.1%, along 
the (001) direction^. Ti02 is a key material for most 
applications, including catalytic and optical devices, sen- 
sors, optoelectronics and spintronics. Anatase thin films 
are frequently grown over LAO by pulsed laser deposi- 
tion. It seems therefore important to perform a careful 
characterization, both from theory and experiments, of 
the film/substrate interface. Just to mention one exam- 
ple where this study may be relevant, room tempera- 
ture ferromagnetism has been obtained from doped and 
undoped anatase films grown over LAO and the results 
strongly depend on the growth conditions^. 

There are two possible interfaces for the system 
LA0/Ti02-anatase, namely LaO facing Ti02 and AIO2 
facing Ti02^'^. Due to the ionicity of the component 
oxides, and in particular, to the fact that the layers of 
LAO (AIO2 and LaO) have alternating formal charges 
(-1 and -|-1, respectively), there is an interfacial formal 
excess charge that should be compensated either by the 
presence of terraces with different stacking, or by oxygen 
vacancies or by atomic interdiffusion. 

Experimental work^ has shown that, in the case of 
LA0/Ti02-anatase, terraces with both types of inter- 
faces appear and approximately in the same proportion. 



For the theoretical modeling of this system the fact that 
the Ti02-anatase layers along the (001) direction are 
not strictly planar must be taken into account. As a 
consequence, the anatase interface can be either oxygen- 
terminated and thus neutral (with formal valence -1-4 for 
the interfacial Ti ion), or Ti-terminated. The latter situa- 
tion is equivalent to considering surface oxygen vacancies 
and gives rise to a change in the valence of the interfacial 
Ti ion. 

In this paper we study the electronic and mag- 
netic properties of interfaces for different stackings of 
LA0/Ti02 using density functional theory^ in the lo- 
cal density approximation (LSDA) and also considering 
electronic correlations in the framework of the LSDA+U 
approximation. 



II. METHOD OF CALCULATION 

To simulate the interfaces, we use a periodically re- 
peated slab-geometry with 5 layers of each of the com- 
ponent materials and enough empty space between slabs 
so that they do not interact with each other (see Fig[T]). 
The cell parameters in the plane of the interface are kept 
fixed to those of the experimental bulk LAO (3.79 A)^ 
and the positions of all atoms are allowed to relax in the 
out of plane direction until the forces arc smaller than 
0.05 cV/ A. 

As in our previous work^, we use an odd number of 
layers of LAO to avoid the formation of a large dipole 
moment, but we are aware that with this procedure it is 
difficult to assure whether the interface will be conduct- 
ing or not. The reason lies in the fact that an LSDA 
calculation for a positively (negatively) charged system 
will locate the Fermi level in the conduction (valence) 
band. 

We use the Wien2k code, that is an implementation of 
the full potential linear augmented plane waves method 
(FP-LAPW)^. The calculations are scalar relativistic 
and the parameters used are listed on Ref. H. It is well 
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FIG. 1: Slab structures for three low energy configurations, (a) Interface LaO/Ti02; (b) Interface A102/TiO, with Ti atom in 
the hollow position; (c) Interface A102/TiO, with Ti atom in the bridge position. Oxygen atoms are co-planar with La and Al 
ions and 0.4 Aoff-plane in the Ti layers. 01 and 02 are below and above the Ti-plane, respectively. Oxygen 02 is closer to 
LAO and it is removed in cases (b) and (c). 



known that in the case of oxides there is a band gap 
underestimation when using local exchange correlation 
functionals within density functional methods. For this 
reason, we have performed calculations with both the 
LSDA and the LSDA+U approximationssiS, using U=0.4 
Ry for the Ti d orbitalaiii, as in previous works. The 
comparison between the two procedures will evidence 
whether there is a qualitative difference in the interfacial 
properties due to electron correlations, as for example, a 
change in the valence of the Ti ions. 



III. RESULTS 

The total energies of different possible structures (at 
OK and with coUinear spins) have been evaluated in pre- 
vious works^, and in this paper we present the results 
for the lowest energy ones. For the interface in which 
LaO faces Ti02 (FiglIJa)), the oxygen-termination for 
anatase is preferred. However, when AIO2 faces Ti02, 
the lower energy is for the Ti-tcrmination, thus with oxy- 
gen vacancies^. Two types of Ti-tcrminated interfaces 
are studied: one where the interfacial Ti atom faces the 
hollow site of the Al atoms in the AIO2 layer (Figlljb)) 
and another one where the Ti atom faces a bridge posi- 
tion of the Al sites in the AIO2 plane (FiglTJc)). 

In all the cases considered there is a total formal charge 
of -fl arising in (a) from the LAO slab while in cases (b) 
and (c) it is due to the contribution of both the LAO slab 
and the oxygen vacancy in the anatase slab. We point out 
that there are more possible structures for each type of 
stacking, and that in this work we consider three of the 
lowest energy ones, those that require more energy to 
separate LAO from anatase. LaO/Ti02 lias the lowest 
energy, followed by A102/TiO where oxygen vacancies 
are presenilis (i.e. 02 is removed in FiglTJb) and (c)). 



In the hollow structure one oxygen atom from the AIO2 
layer moves towards the interfacial Ti thus increasing its 
number of neighbors. This relaxation decreases the total 
energy of the system with respect to the bridge structure 
by 0.2 eV in LSDA and 0.4 eV in LSDA+U. 

FigI2] shows the densities of states (DOS) for the three 
structures depicted in FigHJ We present the results for 
the LSDA-I-U approximation as they are similar to the 
LSDA ones and if there are any differences, they are men- 
tioned in the text. On the left side of Fig 121 we present 
the total DOS and the partial contributions from the La 
atoms. On the right side, the DOS projected on the Ti 
atoms are shown in a zoomed-in region close to the Fermi 
energy (Ey). As in the case of pure anatase, the valence 
band is mostly from the oxygen atoms and the conduc- 
tion band from the Ti atoms, with the d orbitals clearly 
separated by symmetry: the t2g levels closer to E/ and 
the eg ones at higher energies. There are no Al states 
close to E/, but oxygen atoms from LAO contribute sig- 
nificantly to the higher energy part of the valence band. 

LaO/Ti02 interface 

In this interface (FiglIJa)), Ef is in the conduction 
band, as expected due to the formal charge in the LAO 
slab. The largest contribution to the occupied states 
close to Ef comes from the interfacial Ti ion and the 
rest of the Ti atoms contribute less as the distance to the 
interface increases, as shown in Figl^^top-right). This is 
true for both LSDA and LSDA-I-U approximations. It is 
interesting to note that similar layer projected DOS were 
obtained for LAO/STOi^ even though in that structure, 
the Ti02 planes alternate with SrO layers. The authors 
also find that the charge distribution is not affected by 
the inclusion of correlation. There is certainly no contri- 
bution from the interfacial La ions to the occupied states, 
as that peak is shifted towards the right when compared 
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FIG. 2: Densities of states for the three slabs of Fig. [T] in the LSDA+U approximation, (a) Top: LaO/Ti02, (b) Middle: 
A102-hollow, (c) Bottom: A102-bridge. The Fermi level (E/) is marked by the dotted vertical line. Positive and negative 
values correspond to majority and minority spins, respectively. Left: total densities of states of the slabs and La levels. Right: 
Partial DOS on the Ti atoms zoomed-in aromid E/. The labels of the Ti projected DOS are the same for the 3 cases. In the 
legend, the number by the atoms indicate the number of planes away from the interface. 
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FIG. 3: LSDA+U charge density of states close to E/ in the plane defined by the Ti-Ti bonds for one unit cell. The occupied 
states in the energy range: [E/-f .5eV, E/] are plotted for the three structures in Fig[l] namely (a) LaO/Ti02, (b) A102-hollow, 
(c) A102-bridge. Charge isolines are spaced 0.01, from 0.01 to 1.00. In each case, the Ti interfacial atom is placed at the center 
of the figure, so as to make them comparable. 
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FIG. 4: Magnetic moment inside the Ti mufhn-tin spheres 
as a function of the penetration in the anatase slab. The 
interface is at layer and the surface is at layer 4. 



with the La atom at the free surface of the slab (see 
FiglUa)). The charge density at the interface is plotted 
in FiglSlJa) and shows clearly the shape of the Ti d or- 
bitals. The main occupation comes from dxz and dzz2-r2- 
A small total magnetic moment is obtained for the unit 
cell, distributed among all the Ti atoms (see FigU]). 

AlOa/TiO-hollow interface 

This interface (FiglT](b)), presents a difference between 
LSDA and LSDA+U approximations. In LSDA+U, the 
majority spin of the interfacial Ti dxy orbital shifts to- 
wards lower energies, so that the system becomes a semi- 



conductor with a small gap of 0.1 eV (see Figl^Jb)) This 
shift is smaller in the LSDA, so that the system is metal- 
lic. In Figl31[b) we show the charge density at the anatase 
side of the interface where the dxy symmetry is evident. 
Both the charge and the magnetic moment at the Ti ions 
are large at the interface and decrease abruptly inside the 
anatase slab (see FigH]). 

A102/TiO-bridge interface 

In this interface (FigHJc)), E/ lies in the conduction 
band, and the La peaks are further away from it, overlap- 
ping the t2g orbitals of the Ti atoms. The same result is 
obtained with LSDA. The interfacial Ti orbitals have dif- 
ferent symmetry from those of the hollow case, being of 
eg character: dx2-y2 and d'iz2-r2, as can be seen in Figl2] 
and Fig[21 In this case, the Cg levels split into two regions: 
one occupied band close to Ej and another one at ener- 
gies above the t2g orbitals, in the conduction band. The 
extra charge and the magnetic moment are both local- 
ized at the interface, near the oxygen vacancy, and their 
magnitude decreases abruptly inside the anatase slab as 
can be seen in FigU) In this case, the interfacial Ti for- 
mal valence and the magnetic moment are smaller than 
in the hollow structure. 

It is interesting to note that in A102/TiO interfaces, 
there is no contribution of the Ti atom at the free surface 
in the occupied region of the conduction band, showing 
that a slab with 5 layers is a good enough approximation 
to account for interface properties. 

There has been some controversy concerning the ap- 
pearance of interfacial ferromagnetism in the LAO/STO 
interfacei^, which may also be the case here. For this rea- 
son we performed calculations with a double size unit cell 
(54 atoms), for the A102/TiO interfaces which present lo- 
calized magnetism at the interface. We studied the mag- 
netic interaction between the interfacial Ti atoms within 
the LSDA+U method, setting them both in parallel and 
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antiparallel configurations. The two solutions exist but 
the antiparahel one is lower in energy, indicating that a 
long-range ferromagnetic interaction will not arise from 
the considered concentration and distribution of oxygen 
vacancies in A102/TiO interfaces. We do not exclude the 
possibility of ferromagnetism in LaO/Ti02, but it would 
be minor and due to fact that all the Ti atoms have local 
magnetic moments, many magnetic configurations should 
be considered for this study. 

To explore the consequences of a lower vacancy concen- 
tration, and also of a formally neutral system, we consid- 
ered the double size unit cell in the hollow configuration 
with only one interfacial oxygen vacancy. In this calcula- 
tion none of the Ti ions resulted Ti-1-3 and no magnetic 
moments appeared. Thus, a charge imbalance is needed 
to obtain Ti+3 ions and is also a necessary (but not suf- 
ficient) condition for the appearance of magnetism. 

IV. CONCLUSIONS 

The principal conclusion of this work is that whatever 
the origin of the charge at the interfaces, its main effect 
is to change the valence of the Ti atoms, either if it is 
due to the layered-structure of LAO or to the presence 
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of oxygen vacancies in anatase. 

In A102/TiO, there are interfacial oxygen vacancies 
and the interfacial Ti atom close to the vacancy position, 
acquires an extra charge. In LaO/Ti02, there are no 
oxygen vacancies as they are not electrostatically favored 
and the extra charge is distributed among all Ti atoms, 
with decreasing value as the distance from the interface 
increases. Our results indicate that the LaO/Ti02 in- 
terface spans through several anatase layers while the 
A102/TiO one is more localized and presents large local 
magnetic moments. 

In experimental samples, co-existence of the studied 
interfaces and possibly others as well, is expected. As 
a consequence, small patches with different magnetic or- 
der might appear—, thus giving rise to sample-dependent 
results. 
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